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ABSTRACT 

This paper presents the first results of comparison of Planck dust maps at 353, 545 and 857 GHz, along with IRAS data at 3000 (100 yum) and 
5000GHz (60/um), with Green Bank Telescope 21 -cm observations of Hi in 14 fields covering more than 800 deg 2 at high Galactic latitude. The 
main goal of this study is to estimate the far-infrared to sub-millimetre emissivity of dust in the diffuse local interstellar medium (ISM) and in the 
intermediate- velocity (IVC) and high-velocity clouds (HVC) of the Galactic halo. Galactic dust emission for fields with average H i column density 
lower than 2 x 10 20 cm -2 is well correlated with 21 -cm emission because in such diffuse areas the hydrogen is predominantly in the neutral atomic 
phase. The residual emission in these fields, once the H i-correlated emission is removed, is consistent with the expected statistical properties of the 
cosmic infrared background fluctuations. The brighter fields in our sample, with an average Hi column density greater than 2 x 10 20 cm -2 , show 
significant excess dust emission compared to the H i column density. Regions of excess lie in organized structures that suggest the presence of 
hydrogen in molecular form, though they are not always correlated with CO emission. In the higher H i column density fields the excess emission 
at 857 GHz is about 40% of that coming from the Hi, but over all the high latitude fields surveyed the molecular mass faction is about 10%. Dust 
emission from IVCs is detected with high significance by this correlation analysis. Its spectral properties are consistent with, compared to the local 
ISM values, significantly hotter dust (T ~ 20 K), lower sub-millimeter dust opacity normalized per H-atom, and a relative abundance of very small 
grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in 
which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper 
limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and 
depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the 
local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular 
clouds, is active much earlier in the cloud condensation and star formation processes. 

Key words. Methods: data analysis - dust - local insterstellar matter - Galaxy: halo - Sub-millimeter: ISM - Infrared: ISM 
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1. Introduction 

Planck 1 (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third-generation space mission to measure the anisotropy of the 

1 Planck (http://www.QS2L.mi/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
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cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31/ to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0.1 K. Planck's sensitivity, angular resolution, and fre- 
quency coverage make it a powerful instrument for Galactic 
and extragalactic astrophysics as well as cosmology. This paper 
presents the first results of the analysis of Planck observations of 
the diffuse interstellar medium (ISM) at high Galactic latitude. 

From the pioneering work of Spitzer and Field (Spitzer 1956; 
Field 1965; Field et al. 1969), observations of the diffuse ISM in- 
cluding intermediate and high-velocity clouds (IVCs and HVCs) 
have been the basis of our understanding of the dynamical inter- 
play between ISM phases and the disk-halo connection in rela- 
tion to star formation. Space-based observations have given us 
spectacular perspectives on the diffuse Galactic infrared emis- 
sion, which highlight the role of dust not only as a tracer of the 
diffuse ISM but also as an agent in its evolution. 

The Infrared Astronomical Satellite (IRAS) revealed the in- 
tricate morphology of infrared cirrus (Low et al. 1984) and 
prompted a wide range of observations. The cirrus is inferred 
to be inhomogeneous turbulent dusty clouds with dense CO- 
emitting gas intermixed with cold (CNM) and warm (WNM) 
neutral atomic gas and also diffuse H2. From imaging by the 
Spitzer Space Telescope, and very recently by the Herschel 
Space Observatory, their structure is now known to extend 
to much smaller angular scales than observable at the IRAS 
resolution (Ingalls et al. 2004; Martin et al. 2010; Miville- 
Deschenes et al. 2010). Observations from IRAS, the Infrared 
Space Observatory (ISO), and Spitzer have also been used to 
characterize changes in the spectral energy distribution (SED) 
from mid- to far-IR wavelengths, which have been interpreted as 
evidence for variations in the abundance of small stochastically- 
heated dust particles. The correlation with H 1 spectroscopic data 
suggests that interstellar turbulence may play a role in changing 
the dust size distribution (Miville-Deschenes et al. 2002a). 

Since the breakthrough discoveries made with the Cosmic 
Background Explorer (COBE), the study of dust and the diffuse 
ISM structure has also become an integral part of the analysis 
of the Cosmic Microwave Background (CMB) and the Cosmic 
Infrared extragalactic Background (CIB). Our ability to model 
the spatial and spectral distribution of the infrared cirrus emis- 
sion could limit our ability to achieve the cosmological goals of 
Planck, as well as of present balloon-borne and ground-based 
CMB experiments. 

Accordingly, the Planck survey was designed to provide an 
unprecedented view of the structure of the diffuse ISM and its 
dust content. Planck extends to sub-millimeter (submm) wave- 
lengths the detailed mapping of the infrared cirrus by the IRAS 
survey. Its sensitivity to faint Galactic cirrus emission is limited 
only by the astrophysical noise associated with the anisotropy 
of the CIB. The Planck survey is a major step forward from 
IRAS for two main reasons. First, by extending the spectral cov- 
erage to submm wavelengths, Planck allows us to probe the full 
SED of thermal emission from the large dust grains that are the 
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Fig. 1. The median 21 -cm spectrum (dashed) and the standard- 
deviation spectrum (solid) of the Nl field; the shaded back- 
grounds show the LSR velocity ranges used to estimate LVC, 
IVC and HVC components (from light to dark grey). 



bulk of the dust mass. Second, the dust temperatures obtained 
via submm SEDs also help us to disentangle the effects of dust 
column density, dust heating and dust emission cross-section on 
the brightness of the dust emission. 

The scientific motivation of this paper is to trace the struc- 
ture of the diffuse ISM, including its elusive diffuse H2 com- 
ponent, H + components, and the evolution of interstellar dust 
grains within the local ISM and the Galactic halo. We analyze 
the Planck data in selected fields which cover the full range of 
hydrogen column densities from high Galactic latitude cirrus, 
observed away from dark molecular clouds such as, e.g., Taurus 
(Planck Collaboration 201 lu). For all of our fields, we have deep 
21 -cm spectroscopic observations obtained with the Green Bank 
Telescope (GBT). Our data analysis makes use of, and explores, 
the dust/gas correlation by spatially correlating Planck and IRAS 
data with Hi observations. More specifically our study extends 
previous work on the diffuse ISM SED carried out with 7° res- 
olution FIRAS (Far Infrared Absolute Spectrophotometer) data 
(Boulanger et al. 1996) or with 5' resolution 100 yum IRAS data 
(Reach et al. 1998). 

The paper is organized as follows. In Sect. 2 we describe 
the 21 -cm data and the construction of the column density map 
for each H 1 component. In Sect. 3 we describe the Planck and 
IRAS data. Section 4 describes the main analysis of the paper: 
the determination of the Hi emissivities from 353 to 5000 GHz 
(60 to 850 yum). Results are presented in Sect. 5 followed by a 
discussion of some implications in Sect. 6. Conclusions wrap up 
the paper in Sect. 7. 

2. 21 -cm data 

2.1. The Green Bank Telescope cirrus survey 

The 21 -cm Hi spectra exploited here were obtained with the 
100-meter Green Bank Telescope (GBT) over the period 2005 
to 2010, as part of a high-latitude survey of 14 fields (for de- 
tails, see Martin et al., in prep). The total area mapped is about 
825 deg 2 . The adopted names, central coordinates, and sizes of 
the 14 GBT fields are given in Table 1. The brighter fields of 
the sample, the ones that cover the range of column densities 
that spans the H i-H 2 transition, are located in the North Celestial 
Loop region, covering the Polaris flare (POL and POLNOR), the 
Ursa Major cirrus (UMA and UMAEAST), the bridge between 
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Table 1. The H i fields. Columns 2-4 give the central Galactic coordinates and size of each field. Columns 5-10 give the average H i 
column density and its uncertainty (see Section A), and the average LSR velocity and HWHM for each H i component. 



the two (SPIDER) and the interior of the loop (SPC). The largest 
field in our sample (NEP) is centered on the north ecliptic pole, 
a region of high coverage for Planck. Two fields were targeted 
for their specific IVC (DRACO and G86). The faintest fields of 
the sample were selected either because the HVC has a major 
contribution to the total hydrogen column density (AG, MC, SP) 
or because they are known CIB targets (Nl and BOOTES). 

The spectra were taken with on-the-fly mapping. The pri- 
mary beam of the GBT at 21 -cm has a full- width half-maximum 
(FWHM) of 9.T, and so the integration time (4 s) and telescope 
scan rate were chosen to sample every 3.5', more finely than the 
Nyquist interval, 3.86'. The beam is only slightly broadened to 
9.4' in the in-scan direction. Scans were made moving the tele- 
scope in one direction (Galactic longitude or Right Ascension), 
with steps of 3.5' in the orthogonal coordinate direction before 
the subsequent reverse scan. 

Data were recorded with the GBT spectrometer by in-band 
frequency switching, yielding spectra with a local standard of 
rest (LSR) velocity coverage -450 < vlsr < +355 kms -1 at a 
resolution of 0.80 kms -1 . Spectra were calibrated, corrected for 
stray radiation, and placed on a brightness temperature (7b) scale 
as described in Blagrave et al. (2010); Boothroyd et al. (2011). 
A third-order polynomial was fit to the emission-free regions of 
the spectra to remove any residual instrumental baseline. The 
spectra were gridded on the equiareal SFL (Sanson-Flamsteed - 
Calabretta & Greisen 2002) projection to produce a data cube. 
Some regions were mapped two or three times. 

With the broad spectral coverage, all Hi components from 
local gas to HVCs are accessible. The total column density A^hi 
ranges from 0.6 x 10 20 cm -2 in the SP field to 10 x 10 20 cm -2 in 
POL and NEP; the average column density per field ranges from 
1.1 x 10 20 cm" 2 in BOOTES to 6.3 x 10 20 cm" 2 in POL. 

2.2. The H i components 

We use channel maps from the 21 -cm GBT spectra to produce 
maps of the Hi column density in different velocity ranges. For 
convenience we have separated the emission into three compo- 
nents for each field: Low- velocity cloud (LVC), IVC, and HVC. 
The selection of the velocity range for each component is based 
on inspection of both the median 21 -cm spectrum and the spec- 



trum made from the standard deviation of each channel map. An 
example of these two spectra for the Nl field is shown in Fig. 1. 
The dashed and solid lines show the median and standard de- 
viation of the brightness temperature as a function of velocity. 
The standard-deviation spectrum, more sensitive to the structure 
within channel maps, is used here to establish the velocity range 
of the three components in cases where velocity components are 
blended in the median spectrum. The three shaded backgrounds 
in Fig. 1 show the velocity ranges used to calculate the H i col- 
umn density of the three components in this field. 

The brightness temperature of each velocity channel was 
converted to column density, assuming an opacity correction for 
Hi gas with a constant spin temperature T s , to provide an esti- 
mate of the total H i column density of each component: 

N m (x, y) = A x T s £ - In (l - Av, (1) 

where A = 1.823 x 10 18 cm" 2 (K kms" 1 )" 1 . In the optically- 
thin case (assuming T s » 7b) this reduces to Nm(x,y) = 
A2 v 7b(x,v,v)Av. 

We used T s = 80 K which is compatible with the collisional 
temperature found from the H2 observations for column densi- 
ties near 10 20 cm" 2 (Gillmon et al. 2006; Wakker 2006). It is 
also similar to the average H 1 spin temperature (column density 
weighted) found by Heiles & Troland (2003). T s will be higher 
for the WNM, but for these high latitude diffuse lines of sight 
7t> <^c 80 K for the broad WNM lines, and so adopting the wrong 
r s is of no consequence. For most fields the correction is less 
than 3% compared to the optically-thin assumption. Indeed very 
few of our 21 -cm spectra reach brightness temperatures above 
40 K - only 3% of the spectra in POL field, the brightest one in 
the sample. For these "extreme" cases the opacity correction to 
the column density reaches 35%. Figures 2 and 3 show the Hi 
column density maps of all fields, in units of 10 20 cm" 2 . 

2.3. Uncertainties in N m 

The main analysis presented here relies on a correlation analysis 
between far-infrared/submm brightness and A/hi of the compo- 
nents deduced from 21 -cm emission. In order to estimate prop- 
erly the uncertainties of the deduced correlation coefficients, we 
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need to evaluate the uncertainty of the values of TVhi for the H i 
components. The method described in Appendix A takes advan- 
tage of the fact that the GBT observations were obtained in two 
polarisations. The difference between these spectra gives a direct 
estimate of the uncertainty in each channel, which can be inte- 
grated over the appropriate velocity ranges, providing a column 
density uncertainty for each H i component. The average H i col- 
umn densities and uncertainties, expressed in units of 10 20 cm -2 , 
are given in Table 1 for the three H i components of each field. 
This table gives also the average velocity of each H i component 
and an estimate of the half- width at half-maximum of the 21 -cm 
feature. 



3. Planck and IRAS 

3.1. Map construction 

Our analysis uses infrared to submm data at 3000 and 5000 GHz 
(100 and 60 /mi, respectively) from IRAS (IRIS, Miville- 
Deschenes & Lagache 2005) and at 353, 545, and 857 GHz 
(850, 550, and 350/mi, respectively) from Planck (DR2 release; 
Planck HFI Core Team 201 lb), beginning with maps in Healpix 
form (Gorski et al. 2005) with Af side = 2048 (pixel size of l.V). 
We concentrated here on the three highest frequencies of Planck 
to avoid the significant contamination from residual CMB fluc- 
tuations and interstellar emission other than thermal dust (CO, 
synchrotron, free-free and spinning dust). 

To obtain infrared- submm maps corresponding to each GBT 
field we first projected each Healpix map, using the nearest 
neighbour method, onto SFL grids with a pixel size of 1.7'. Each 
grid was centred on a given GBT field with a size 10% larger in 
each direction in order to avoid edge effects in subsequent con- 
volution steps. Each SFL map was then converted to MJysr -1 
and point sources were removed and replaced by interpolation 
of the surrounding map 2 . The map was then convolved to bring 
it to the GBT 9.4' resolution and finally projected, using bi-linear 
interpolation, on the actual GBT grid (3.5' pixel -1 ). The Planck 
and IRAS maps for our fields are shown in Figs 4 to 1 1 . As will 
be discussed in Sect. 4, these figures also show models of this 
emission based on Hi observations in a masked subset of the 
map, and the residual map on subtracting this model from the 
entire field. The residuals are largest for those areas in the map 
not used to constrain the model. 

3.2. Dust brightness uncertainty 

To estimate the noise level of the Planck and IRAS maps we 
used the method described in Sect. 5.1 of Miville-Deschenes & 
Lagache (2005). For both data sets we used the difference of 
maps of the same region of the sky obtained with different sub- 
samples of the data. These difference maps, properly weighted 
by their coverage maps, provide an estimate of the statistical 
properties of the noise. For Planck the noise was estimated us- 
ing the difference of the first and second half ring maps (Planck 
HFI Core Team 2011b). In the case of IRAS, each ISSA plate is 
the combination of up to three maps built from independent ob- 
servations over the life of the satellite. We built difference maps 
from these three sets of maps. The procedure used to estimate 
the Planck and IRAS noise levels at the GBT resolution is de- 



2 For Planck channels we removed only point sources identified in 
the ERCSC (Planck Collaboration 2011c). For IRAS maps we used 
the source removal method described in Miville-Deschenes & Lagache 
(2005). 



tailed in Appendix B. The noise levels for each field and each 
frequency are given in Table B.l. 

4. Dust-HI correlation 

4.1. Model 

Many studies, mostly using the IRAS and COBE data com- 
pared with various 21 -cm surveys (Boulanger & Perault 1988; 
Joncas et al. 1992; Jones et al. 1995; Boulanger et al. 1996; 
Arendt et al. 1998; Reach et al. 1998; Lockman & Condon 2005; 
Miville-Deschenes et al. 2005), have revealed the strong corre- 
lation between far-infrared/submm dust emission and 21 -cm in- 
tegrated emission Wm 3 at high Galactic latitudes. In particular 
Boulanger et al. (1996) studied this relation over the whole high 
Galactic latitude sky. They reported a tight dust-H i correlation 
for Wm < 250 K km s" 1 , corresponding to N m < 4.6x 10 20 cm" 2 . 
For higher column densities the dust emission systematically ex- 
ceeds that expected by extrapolating the correlation. Examining 
specific high Galactic latitude regions, Arendt et al. (1998) and 
Reach et al. (1998) found infrared excesses with respect to TVhi, 
with a threshold varying from 1.5 to 5.0 x 10 20 cm . 

Part of this excess is due to the effect of 21 -cm self- 
absorption that produces a systematic underestimate of the col- 
umn density when deduced with the optically thin assumption. 
Even though this effect is only at a level of a few percent in our 
case because of the low column densities, applying an opacity 
correction (see Sect. 2.2) helps to limit this systematic effect. 

Most of the infrared/submm excess is usually attributed to 
dust associated with hydrogen in molecular form. This hypothe- 
sis is in accordance with UV absorption measurements that show 
a sudden increase of the H2 absorption at Nu = (3-5)xl0 20 cm -2 
(Savage et al. 1977; Gillmon et al. 2006), roughly the threshold 
for departure from the linear correlation between dust emission 
and A^hi- It is also observed that the pixels showing evidence of 
excess are spatially correlated and correspond to, or at least are 
in the vicinity of, known molecular clouds traced by CO emis- 
sion. See also the discussion in Sect. 6.1. 

A third source of this excess emission could be dust asso- 
ciated with the Warm Ionized Medium (WIM) but detection of 
this component is difficult (Arendt et al. 1998; Lagache et al. 
2000) because there is no direct tracer of the ionized gas column 
density; Ha depends on the square of the electron density and 
part of the structure seen in Ha might be back- scattering of dif- 
fuse Galactic emission on dust and not photons produced within 
cirrus clouds (Witt et al. 2010). 

Finally, in the most diffuse regions of the high-latitude sky, 
the fluctuations of the CIB are a significant fraction of the bright- 
ness fluctuations in the infrared/submm. With a power spec- 
trum flatter (k~ l ) than that of the interstellar dust emission (k~ 3 ) 
(Miville-Deschenes et al. 2002b, 2007; Lagache et al. 2007; 
Planck Collaboration 201 In), the CIB anisotropies contribute 
mostly at small angular scales, producing statistically homoge- 
neous brightness fluctuations over any observed field, like an in- 
strumental noise. Furthermore, because the CIB is unrelated to 
interstellar emission, the CIB fluctuations cannot be responsible 
for the excess of infrared emission seen at moderate to high Nui 
column density. 

In the analysis presented here we go a few steps further than 
the previous studies by: 1) allowing for different dust emissivi- 
ties for the local ISM (i.e., LVC), IVC, and HVC components; 
2) applying an opacity correction to the 21 -cm brightness tem- 
peratures in order to compute a more reliable A/hi; and 3) taking 

3 equivalent to optically-thin H 1 column density. 
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Fig. 6. Like Fig. 4, for field SPC. 



into account explicitly the CIB fluctuations which turn out to 
dominate the uncertainties in the derived emissivities. 

In some fields (like G86 with strong IVC emission) the dis- 
tinctive morphology of the IVC column density map can be seen 
clearly in the line-of-sight integrated dust emission map (see 
Figs. 2 and 5 and Martin et al. 1994), but even faint signals can 
be brought out by formal correlation analysis. We use the fol- 
lowing model: 



where I v (x, y) is the dust map at frequency v (IRAS or Planck), e l v 
the emissivity of Hi component i (LVC, IVC and HVC), and Z v 
is the zero level of the map. R v represents not only the contribu- 
tion from noise in the data but also any emission in the IRAS and 
Planck bands that is not correlated with Nui including the CIB 
anisotropics and potential dust emission coming from molecular 
or ionized gas. In this model we assume that the three HI com- 
ponents N l m (x,y) have a constant emissivity e\ over the field. 
Any spatial variations of the emissivity would also contribute to 
fluctuations in R v (x,y). 



(2) 
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Fig. 7. Like Fig. 4, for field UMAEAST. 



4.2. Estimating the dust emissivities 

To estimate the parameters e\ and constant Z v we used the IDL 
function regress which, in the case of a general linear least- 
squares fit, solves the following equation (Press et al. 1995): 



a = (A T A)~ l x (A T b), 



(3) 



where a is the vector of the parameters e l v and b is a vector of 
the N IRAS or Planck data points from the map, divided by their 
respective error: 



(4) 



A is an N x M matrix that includes the Nui values of the M H i 
components, 



*J = 



(5) 



Regress uses a Gaussian elimination method for the inversion. 

For the model described in Eq. 2, the least-squares fit method 
provides a maximum-likelihood estimation of the parameters e\ 
provided that the residual term R v (x, y) is uncorrected with N l m 
and its fluctuations are normally distributed (i.e., white noise). 



In addition, in order that the parameter estimates not be biased, 
the uncertainties on N l m have to be comparatively small; we will 
show (see Sect. 4.5) that this last condition is satisfied for our 
data. On the other hand, we will also show that the residual 
term R v (x, y) is clearly not compatible with white noise. Even 
for the most diffuse fields in our sample, where the CIB fluctua- 
tions dominate the residual emission and the Probability Density 
Function (PDF) of R v (x, y) is normally distributed, the condition 
that there be no (not even chance) correlation with TVhi is not sat- 
isfied, because the power spectrum of the CIB is not white, but 
rather like k~ l (Planck Collaboration 201 In). For brighter fields, 
where spatial variation of the dust emission with respect to the 
Hi templates is expected (due to the presence of molecules, a 
poor H i opacity correction, or spatial variation of dust proper- 
ties), the residual is not even normally distributed. 



To limit the influence of these effects, and in order to focus 
on estimating the dust emissivity of the H i components, we re- 
lied on a masking procedure to flag and remove obvious outliers 
with respect to the correlation, and on Monte-Carlo simulations 
to estimate the e l y uncertainties and bias (see Sect. 4.5). 
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Fig. 8. Like Fig. 4, for field UMA. 



4.3. Masking 

In order to limit the effect of lines of sight with significant "ex- 
cess" dust emission that is not associated with H i gas, most pre- 
vious authors used only data points with Nm lower than a given 
threshold to stay in a regime of linear correlation. This thresh- 
olding was motivated by the fact that above some Nm the extinc- 
tion and self- shielding of H2 are strong enough to limit photo- 
dissociation, whereas below the threshold the hydrogen is mostly 
atomic. The threshold depends sensitively on local gas density 
and temperature but for typical interstellar conditions for CNM 



gas (n = 100 cm" 3 , T = 80 K, G = 1, where G is the scaling 
factor of the Interstellar Radiation Field (ISRF) as defined by 
Mathis et al. (1983)), it is about Af ffl = 2.5 x 10 20 cnT 2 (Reach 
et al. 1998). Others have used a quadratic function for e y (Fixsen 
et al. 1998) based on the idea that the H2 column density depends 
(at least dimensionally) on (Reach et al. 1994). Both meth- 
ods introduce a bias in the parameter estimation that is difficult 
to quantify. 

Instead of applying an arbitrary cutoff in Nm, Arendt et al. 
(1998) used an iterative method to exclude data points above a 
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cut along lines perpendicular to the fit, in order to arrive at a correspond to lines of sight where the dust emission is dominated 
stable solution for e v . We used a similar approach by iteratively by the atomic Hi components, 
masking out data points that would produce a positively- skewed 
residual. That way we expect to keep pixels in the maps that 
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The PDF of residual map, R, defined as 



R v (x,y) = I v (x,y) - ^ e v N l m (x,y) - Z v 
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(6) 



is used to estimate the mask. We used the Planck 857-GHz chan- 
nel because it has the best signal-to-noise ratio and is less sen- 
sitive than the IRAS channels to dust temperature-induced emis- 
sivity variations. For the most diffuse fields in our sample (AG, 
MC, Nl, BOOTES, G86, SP) the PDF of R S51 is very close to a 



13 



Planck collaboration: Planck Early Results. XXIV. Dust in the diffuse interstellar medium and the Galactic halo 



NEP 




100° 96° 92° 100° 96° 92° 100° 96° 92° 

Gal. Ion. Gal. Ion. Gal. Ion. 



11. Like Fig. 4, for field NEP, the largest. 



Planck collaboration: Planck Early Results. XXIV. Dust in the diffuse interstellar medium and the Galactic halo 



1.0 

0.5 

0.0 
1.0 

0.5 

0.0 
1.0 

0.5 

0.0 
1.0 

0.5 

0.0 
1.0 

0.5 

0.0 
1.0 

0.5 

0.0 
1.0 

0.5 





\ BOOTES 
\ a=-0.70 




■ / 


\ N1 
I a=0.07 




■ / 


1 G86 
\ a=0.09 






v DRACO 
1 a=0.08 



SP 
a=0.07 



AG 
a=0.08 





\ MC 


/ 


1 a=0.07 




SPC 
a=0.10 




NEP 
a=0.12 




0.0 







UMAEAST 






a=0.12 











POL 




a=0.15 







-1.0 -0.5 0.0 0.5 1.0 1.5 

R 857 (MJysr 1 ) 



-1.0 -0.5 0.0 0.5 1.0 1.5 
R 857 (MJy sr" 1 ) 



Fig. 12. Normalised PDF of the residual R^sj of the dust-gas 
correlation at 857 GHz for each field after convergence of the 
masking procedure. The red curve is the result of a Gaussian fit 
to the lower, rising part of the PDF. 



Gaussian, which suggests that the model described by Eq. 2 is 
the right one in such low column density regions. 

Accordingly, for the first iteration of the masking process 
for each field, we performed the multi-variate linear regression 
based on Equations 3 to 5 using only the faintest 10% pixels 
in the 857 GHz map. For further iterations this threshold was 
relaxed, bringing in more pixels in the map compatible with the 
iterated model. 

As discussed above, the presence of dust emission associated 
with molecular gas can positively skew the PDF, and empirically 
the PDF of 7?857 is indeed positively skewed for the eight remain- 
ing fields (Fig. 12). To determine the set of pixels to be retained, 
we used a Gaussian fit to the lower, rising part of the PDF, up to 
the PDF maximum, and estimated the cr (see the red curves in 
Fig. 12). With the above motivation, we assume that the lower 
part of the PDF is representative of pixels where the fit works 
well (i.e., for these pixels the residual is normally distributed). 
Using the cr fit only to this part of the PDF, we apply a threshold 
in 7?857 by masking out all pixels with 7? 8 5 7 > 3 x cr away from 
the mean. We iteratively recompute the parameters and mask to 
converge on a stable solution. The PDFs of R^ obtained at the 
end of the process are shown in Fig. 12. 

For the six faintest fields, the masking excluded less than 1% 
of the points. For these fields the mask has no significant ef- 
fect on the estimated parameters. For the eight other fields, the 
masking method excluded from 17 to 83% of the pixels. The 
masks for these eight fields are shown in Fig. 13. In these cases 
the masking has a significant effect on the result, but we have 



checked that the estimated parameters are similar to the ones ob- 
tained with a A^hi < 4 x 10 20 cm -2 threshold. In fact the masking 
method used here allows us to keep pixels that would have been 
excluded by a simple Nui thresholding even though they do not 
depart significantly from the linear correlation. 

Table 2 provides the e l v values for each 
field/component/frequency. In order to visualize the re- 
sults, Fig. 14 gives scatter plots together with the line of slope 
6g 5? for each component, field by field. Specifically, for each HI 

component i we plotted I v - e jN J m - Z y as a function of N l m . 



4.4. Statistics of the residua! 

Figs 4 to 11 show the IRAS and Planck maps, together with 
the H i correlated emission and the residual maps R v for all our 
fields. For the six faintest fields (Nl, SP, BOOTES, AG, MC, and 
G86 - see Figs. 4-5), the structure in the residuals is, even by vi- 
sual inspection, clearly spatially correlated between frequencies, 
especially in the Planck bands. It is dominated by small scale 
structures with equally negative and positive brightness fluctu- 
ations. The structure of the residual for brighter fields is also 
clearly correlated between frequencies but in these cases the 
residual is mostly positive (i.e., they are excesses with respect 
to the Hi). These residuals also show larger coherent structures 
than in the fainter fields. 

The rms (about the mean) of the residual R can be approxi- 
mated as: 



ctr 



+ (cr dust ) 2 + 

v noise 7 



(7) 



where noise in the IRAS or Planck data and that induced by un- 
certainties in A^hi are explicitly accounted for and cr s includes all 
other contributions not in the model, including CIB anisotropies 
and dust emission associated with molecular gas. After quadratic 
subtraction Fig. 15 shows the value of cr s at each frequency, as a 
function of the average Nm density for each field. For fields with 
a median column density lower than 2 x 10 20 cm~ 2 , the PDFs 
of the residual emission all have skewness and kurtosis values 
compatible with a Gaussian distribution. Furthermore the width 
of these PDFs shows very small scatter from field to field (see 
Fig. 15). This is another indication that, for such diffuse fields, 
the model is a good description of the data; i.e., the Galactic dust 
emission is largely dominated by the H i components, with very 
limited spatial variations of the emissivity across a given field. 

The dashed line gives the average level of cr s for those fields 
with TVhi < 2 x 10 20 cm~ 2 . These quantities, together with the 
average Planck or IRAS and GBT noise contributions to cr R , are 
summarized in Table 3. It is clear from Table 3 and Fig. 15 that 
even for faint fields cr s is the main contributor to the rms of the 
residual, and therefore the main contributor to the dispersion in 
the dust-gas correlation diagrams in these fields. This is in ac- 
cordance with the findings of Planck Collaboration (201 In) who 
concluded that the CIB fluctuations dominate os at 353 GHz and 
higher frequencies. 

At higher column densities, the PDF of the residual shows 
positive skewness and an rms that increases with Nm, signifi- 
cantly exceeding the level of CIB anisotropies. These aspects of 
the PDF indicate one or more extra components that contribute 
to the dust emission and are not taken into account in our model. 
Contributions to the residual that grow with Nui are compatible 
with the presence of a molecular gas component. This could also 
come from spatial variations of the dust emissivity for given H i 
components or from inadequately-separated Hi components. 
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Fig. 13. Mask for the eight fields where more than 1% of the pixels were excluded (NEP=17%, POL=83%, SPC=26%, 
POLNOR=30%, SPIDER=24%, UMA=56%, UMAEAST=66%, DRACO=29%). For each field the left image is the total Hi 
integrated emission and the right image is the Planck 857 GHz map. The regions in greyscale were excluded from the correlation 
analysis. Note the relationship of the masks to the regions of high residual in Figs. 5 to 11. 



16 



Planck collaboration: Planck Early Results. XXIV. Dust in the diffuse interstellar medium and the Galactic halo 



Field 


HI 


^353 


^545 


^857 


£3000 


^5000 


AG 


LVC 


A no A r\ A A -"7 

0.034 + 0.007 


M 1 A A AO 

0.14 + 0.02 


A O M , M MC 

0.39 + 0.05 


A f Ci M M A 

0.68 + 0.04 


A 10 1 , A M 1 O 

0.181 + 0.018 




IVC 


m mom . m r\f\c 

0.020 + 0.005 


a Anr i a m 1 

0.075 + 0.016 


A O O A AO 

0.22 + 0.03 


M C O . M MO 

0.58 + 0.03 


A 1/^1 , MMIO 

0.161 + 0.013 




HVC 


M MM A i M MM1 

0.004 + 0.003 


A M1 1 i M f\C\0 

0.011 + 0.008 


MMIO i MM1C 

0.018 + 0.015 


M MO n i M M1 C 

0.037 + 0.015 


M MMM 1 M AA/T 

0.009 + 0.006 


BOOTES 


LVC 


f\ f\ A C . A A1 A 

0.045 + 0.010 


0.17 + 0.03 


M C A i M M f 

0.54 + 0.06 


M M C . A MO 

0.95 + 0.08 


MOO 1 M MO 

0.23 + 0.03 




IVC 


0.029 + 0.007 


A 1 O . A AO 

0.12 + 0.02 


M T7 i MM/1 

0.37 + 0.04 


M O n . m MC 

0.87 + 0.05 


M T7C , M MOM 

0.275 + 0.020 




HVC 


A AO r\ AO 

0.02 + 0.02 


A AA . A AH 

0.09 + 0.07 


M OM M 1 O 

0.20 + 0.13 


MOM 1 M 1 O 

-0.30 + 0.18 


M 1 O . M Mn 

-0.13 + 0.07 


TNT) A /"< i~\ 

DRACO 


LVC 


A A/jO . M M1 1 

0.043 + 0.011 


A 1 o i M M A 

0.18 + 0.04 


M A C\ i M MT 

0.49 + 0.07 


m cn 1 m Mn 

0.57 + 0.07 


M MO 1 M MO 

0.08 + 0.02 




IVC 


A A /I O r\ AM A 

0.042 + 0.004 


A 1 z' O , A A 1 O 

0.168 + 0.013 


A A O A AO 

0.48 + 0.03 


A nM A AO 

0.70 + 0.03 


A 1 S n A AAA 

0.167 + 0.009 




HVC 


A r\r\i-i , r\ AM f 

0.007 + 0.006 


A AO O , A f\ 1 O 

0.032 + 0.018 


A A^7 , M M A 

0.07 + 0.04 


M 1 M 1 M M A 

0.10 + 0.04 


M MO C ■ M M1 O 

0.025 + 0.013 


G86 


LVC 


a ao o m mm a 

0.033 + 0.004 


A \ A Z' A A 1 O 

0.146 + 0.012 


A AC A AO 

0.45 + 0.02 


0.71 + 0.03 


A 1 A C , A MMO 

0.165 + 0.008 




IVC 


m m 1 c 1 . m mm 1 m 

0.0151 + 0.0019 


0.070 + 0.006 


MOOO i MMIO 

0.238 + 0.013 


M/^/lO 1 A M 1 C 

0.643 + 0.015 


A O A , A AM A 

0.206 + 0.004 




HVC 


m f\r m mo 

-0.05 + 0.02 


A 1 A A AH 

-0.19 + 0.07 


-0.36 + 0.16 


M O A M 1 M 

-0.24 + 0.19 


A A A A MC 

-0.09 + 0.05 


MC 


LVC 


m mo 1 . m m 1 m 

0.031 + 0.010 


A 1 C , A MO 

0.15 + 0.03 


M A C i M M/^ 

0.45 + 0.06 


a n . a An 

0.78 + 0.07 


M 1 O 1 M MO 

0.18 + 0.03 




IVC 


A AMD M MMM 

0.008 + 0.009 


A AO A AO 

0.02 + 0.03 


A 1 A MC 

0.16 + 0.05 


A nM A A f 

0.70 + 0.06 


A 1 n A AO 

0.17 + 0.03 




HVC 


m r\r\s~ . m mmo 

-0.006 + 0.002 


a ao m . a r\r\ r i 

-0.020 + 0.007 


MM/IT i MM1C 

-0.047 + 0.015 


A MOM ■ M M 1 n 

-0.030 + 0.017 


A M 1 M . A MMn 

0.019 + 0.007 


Nl 


LVC 


A nr/ r\ r\f\i—t 

0.056 + 0.007 


A <-| 1 f A AO M 

0.215 + 0.020 


A C O ■ M M A 

0.58 + 0.04 


M O S A AO 

0.86 + 0.03 


A \ f f AMI 1 

0.166 + 0.011 




IVC 


A AO r\ , f\ AAH 

0.039 + 0.007 


A 1 C , A MO 

0.15 + 0.02 


M A 1 i M M A 

0.41 + 0.04 


A n O A AO 

0.72 + 0.03 


moio i mmio 
0.213 + 0.012 




HVC 


M MMC M AA A 

0.005 + 0.004 


A A 1 C A A 1 /I 

0.015 + 0.014 


M M A . M AO 

0.04 + 0.03 


AMI A AO 

-0.01 + 0.02 


A AM 1 A A An 

-0.001 + 0.007 


NEP 


LVC 


M M/IOM . M MM 1 A 

0.0420 + 0.0014 


A 1 A O , A f\f\C 

0.163 + 0.005 


M A ^1 f\ , M M1 1 

0.470 + O.Oll 


A f f A , MMIO 

0.664 + 0.013 


M 1 A 1 . M MMC 

0.141 + 0.005 




IVC 


A A1 C\ r l M AA 1 O 

0.0197 + 0.0012 


A A O A A A A A 

0.080 + 0.004 


A O O f , A A 1 A 

0.236 + 0.010 


A /'/'/' A A 1 O 

0.666 + 0.013 


A O O A A A A A 

0.229 + 0.004 




HVC 


m mo i a mmm 

-0.021 + 0.009 


A AM , M MO 

-0.09 + 0.03 


A O O A f\ r l 

-0.22 + 0.07 


M A O 1 M 1 M 

-0.43 + 0.10 


A AO A AO 

-0.02 + 0.03 


POL 


LVC 


A AC 1 A A A A 1 A 

0.0519 + 0.0019 


A OMO A A AT 

0.203 + 0.007 


A n A AO 

0.57 + 0.02 


A A C C MMIO 

0.455 + 0.018 


A 1 AO A A A A 

0.102 + 0.004 




IVC 


A (\C f , A A1 O 

0.056 + 0.012 


M O A i M MC 

0.24 + 0.05 


0.61 + 0.13 


m no M 1 1 

0.72 + 0.11 


M 1 O 1 M MO 

0.18 + 0.03 




HVC 












POLNOR 


LVC 


m m a 'i s~ . a mm 1 o 

0.0476 + 0.0012 


A O MM . A AM A 

0.200 + 0.004 


A S 1 O , A M 1 O 

0.612 + 0.012 


MCOO ■ A M 1 O 

0.538 + 0.012 


A A O O . A AM O 

0.088 + 0.003 




IVC 


A MOO A C\f\ r l 

0.023 + 0.007 


A MO A AO 

0.08 + 0.02 


A O M A AH 

0.20 + 0.07 


A C 1 A An 

0.51 + 0.07 


A 1 C S . A A 1 f 

0.156 + 0.016 




HVC 












SP 


LVC 


A A A O A A A O 

0.063 + 0.008 


A O C A AO 

0.25 + 0.03 


A no A MC 

0.72 + 0.05 


M CM M M A 

0.59 + 0.04 


A A A A , A A 1 f 

0.094 + 0.016 




IVC 


a AiA , a mm o 

0.029 + 0.008 


M 1 1 ■ M MO 

0.11 + 0.02 


A OM A MC 

0.29 + 0.05 


m nn 1 M M A 

0.77 + 0.04 


A O O M M M 1 f 

\).229 + 0.016 




HVC 


A AAO A A AT 

-0.003 + 0.007 


A AO A AO 

-0.02 + 0.02 


A An M M A 

-0.07 + 0.04 


M 1 O M M A 

-0.12 + 0.04 


A A /I O A M 1 O 

-0.042 + 0.013 


SPC 


LVC 


A f C , A MM 1 i~l 

0.0365 + 0.0017 


M 1 A M i M MM C 

0.140 + 0.005 


A /I M 1 , MM11 

0.401 + 0.011 


A /111 , MM11 

0.4H + 0.011 


A A O A . A AM A 

0.086 + 0.004 




IVC 


A A 1 C A A A O 

0.015 + 0.008 


A A A , A AO 

0.06 + 0.02 


A 1 O A MC 

0.18 + 0.05 


M C A A MC 

0.54 + 0.05 


A 1 An A AO M 

0.197 + 0.020 




HVC 


a mm a . a mm c 

-0.004 + 0.005 


A AM C , M M 1 T 

-0.005 + 0.017 


A AM , M M A 

0.00 + 0.04 


M MC i M M A 

-0.05 + 0.04 


M MMC 1 M M1 O 

-0.005 + 0.013 


SPIDER 


LVC 


M M A <"7 A A A A 1 /I 

0.0474 + 0.0014 


A O MA A A A A 

0.200 + 0.004 


A A A O A A 1 O 

0.602 + 0.012 


A CnM A A 1 O 

0.570 + 0.013 


A A A O A A A A 

0.093 + 0.004 




IVC 


A MOM . A MM O 

0.030 + 0.003 


A 1 f\>~l , A M 1 O 

0.107 + 0.012 


MOM i M MO 

0.30 + 0.03 


M CM . M MO 

0.59 + 0.03 


A 1 f O , A AMM 

0.162 + 0.009 




HVC 


A f\C f A A 1 

-0.056 + 0.016 


M 1 T A AC 

-0.17 + 0.05 


M CO M 1 A 

-0.52 + 0.14 


A OM A 1 

-0.80 + 0.16 


M 1 M M M A 

-0.10 + 0.04 


UMA 


LVC 


0.049 + 0.002 


0.2H +0.007 


0.62 + 0.02 


0.563 + 0.017 


0.098 + 0.005 




IVC 


0.030 + 0.005 


0.12 + 0.02 


0.37 + 0.06 


0.58 + 0.05 


0.151 +0.012 




HVC 


-0.013 + 0.005 


-0.064 + 0.019 


-0.18 + 0.05 


-0.12 + 0.04 


-0.006 + 0.011 


UMAEAST 


LVC 


0.031 +0.003 


0.147 + 0.008 


0.48 + 0.02 


0.566 + 0.016 


0.106 + 0.005 




IVC 


0.059 + 0.004 


0.219 + 0.014 


0.59 + 0.04 


0.65 + 0.03 


0.191 +0.009 




HVC 


0.007 + 0.004 


0.024 + 0.013 


0.06 + 0.03 


0.04 + 0.03 


0.008 + 0.008 



Table 2. Emissivities of each H 1 component at 353, 545, 857, 3000, and 5000 GHz. Units are MJy sr VlO 20 cm 2 . The uncertainties 
were obtained using Monte-Carlo simulations (see Section 4.5). 
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353 


0.0120 + 0.0006 


0.0060 + 0.0012 


0.0014 + 0.0006 


545 


0.038 + 0.002 


0.0096 + 0.0018 


0.006 + 0.003 


857 
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0.0097 + 0.0019 


0.016 + 0.007 


3000 
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0.0146 + 0.0016 


0.0053 + 0.0010 



Table 3. Average of the standard deviations of the sky residual 
cr s for the six fields with (A^hi) lower than 2 x 10 20 cm" 2 (AG, 
BOOTES, G86, MC, Nl, and SP). cr dust is the average level of 

' ' ' ' J noise b 

noise rms in the IRAS and Planck maps at the GBT resolution 
(i.e., noise level in maps convolved at 9.4 r ) computed for the 14 
fields, cr^oise gi yes me average GBT noise converted to MJy sr -1 
(see Eq. 7), for the 14 fields. All uncertainties are the Icr value 
of each sample. 



4.5. Monte-Carlo method to determine emissivity 
uncertainties 

The statistical uncertainties estimated for the emissivities by the 
least-squares fit method are accurate only for the case of white 
Gaussian noise in I v and (sufficiently) low noise in the A^hi com- 
ponents. The noise in I v includes IRAS or Planck instrumental 
noise, CIB anisotropics, and various interstellar contaminants, 
and for the most diffuse fields in our sample, its PDF is close to 
Gaussian. However, its power spectrum is certainly not white. 
First, at the angular scales of our observations (from 9 r to a 
few degrees), the power spectrum of the CIB anisotropics is 
P(k) oc k~ l (Planck Collaboration 201 In). Second, the spatial 
variation of the coverage and the convolution to the GBT reso- 
lution both introduce spatial structure in the noise that modifies 
its power spectrum. The addition of all those noise sources pro- 
duces a net noise term n v on I v that is not white. Because of the 
random chance correlation of n v with the A^hi components, the 
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Fig. 14. 7 8 57 vs. A^hi scatter plots, visualising the dust-gas correlation for each of the three Hi components, across a row. The 
fields POL and POLNOR do not have an HVC component. For a given H i component, the remaining 857 GHz emission, once 
the contribution of the other two Hi components has been removed, is plotted as a function of A/hi of that component (i.e., I v - 
Yij^i ~ Z v as a function of N l m ). The contours show quartiles of the density of points in each scatter plot. The filled contours 
show the data points used in the correlation analysis (i.e., data points not masked out) while the open contours show all data points. 



uncertainties on the parameters estimated from the least- squares 
are systematically underestimated (the least-squares fit is not op- 
timal). In addition, the Nui maps are not noise free; significant 
noise on the independent variable in a least-squares fit produces 
a systematic bias in the solution. Furthermore, imperfect opacity 
correction of the H i spectra, the presence of molecular gas, and 
spatial variations of the dust properties will also produce non- 
random fluctuations in the residual map. For all of these reasons, 



an analysis of Monte-Carlo simulations is required for proper 
estimation of the uncertainties and biases in the e\. 

To generate simulations for each field, we adopted the Nm 
maps obtained from the 21 -cm observations as templates of the 
dust emission. We built dust maps V y for each frequency v by 
adding up these A^hi maps multiplied by their respective es- 
timated emissivities e\ (just as in computing the residual), to 
which we added realizations of the IRAS or Planck noise n v ap- 
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Fig. 15. Value of crs obtained from the standard deviation cr R 
of the residual maps R v (x, y) from which contributions from the 
IRAS or Planck and the GBT noise were removed quadratically 
(see Eq. 7), plotted as a function of the average HI column den- 
sity of each field (sum of LVC, IVC and HVC). The dashed line 
is the average of cr s for fields with A/hi < 2 x 10 20 cm -2 (see 
Table 3). 



propriate to the field 4 at a level compatible with Table 3 once 
convolved at GBT resolution, and of the sky residual a v with a 
k~ l power spectrum at a level to reproduce cr s within the mask 5 : 



some (undetected) HVCs, the bias is only at the few percent 
level; in what follows we made no correction for it. 



5. Results 

5.1. Planck and IRAS emissivities 

The present study extends to smaller scales, and to the IVCs and 
HVCs, the earlier work done with the FIR AS (Boulanger et al. 
1996) or IRAS (Boulanger & Perault 1988; Reach et al. 1998) 
data on the dust emission of the diffuse ISM. It also extends 
to a much larger sample a similar analysis done on a diffuse 
3° x 3° region at high Galactic latitude using IRAS and Spitzer 
data (Miville-Deschenes et al. 2005). The IR/submm-Hi corre- 
lation analysis allows us to determine empirically the spectral 
dependence of the ratio between the dust emission and the gas 
column density. In addition, the combination of Planck, IRAS, 
and A^hi data can be used to trace one elusive component of the 
diffuse interstellar medium: the diffuse H 2 gas (Sect. 6.1). 

As seen in Fig. 14 there is a clear correlation between the 
IRAS or Planck data and Nui in all fields but, as previously 
seen with COBE and IRAS data, there are increasing excesses 
of dust emission with increasing Nw The estimated emissivities 
for each field/component/frequency are compiled in Table 2 and 
shown in Fig. 16. We note that dust associated with the LVC and 
IVC components is detected in each field and at each frequency, 
unlike for HVCs for which we do not report any significant de- 
tection. All HVC emissivities are indeed below 3<x. The results 
on HVCs are discussed further in section 6.3. In the following 
we analyse what can be drawn from the emissivities for the LVCs 
and IVCs. 



Yj 6 ^hi + a v 



+ n v 



(8) 5.2. Comparison with FIRAS data 



Before fitting these simulated data we added white noise to each 
A^hi component at the level estimated for the GBT data (Table 1). 
We then carried out the least- squares fit, using the mask already 
estimated for each field. From a thousand such simulations for 
each field, we obtained the statistics of the recovered ej,' so that 
we could determine if our original fit (fed into the simulation) 
was biased and could compare the dispersion of the parameters 
to the statistical uncertainties returned by the least- squares fitting 
procedure. 

Table 4 summarizes the results of the Monte-Carlo simula- 
tions. Due to random correlation between n v and the A^hi compo- 
nents, we find the Monte-Carlo-derived uncertainty in e l v is sev- 
eral times higher than the analytically-derived uncertainty. The 
results are reported in the table in terms of y , the ratio of <r(6^), 
the standard deviation of the emissivities recovered in the sim- 
ulations, to 5e l v , the standard deviation expected from the linear 
fit performed by Regress. Therefore, in what follows we use the 
values of <x(£") as the uncertainties in e l v (see Table 2). 

Finally, these simulations make it possible for us to estimate 
any bias in e l y which could arise from the noise in A/hi. Table 4 
also provides the bias in %: b v = 100((6^) - e v )/e v . Except for 



4 We assumed white noise maps for both IRAS and Planck, each di- 
vided by the square root of their coverage map. 

5 For the faintest fields this will produce a map with the statistical 
properties of the CIB anisotropies. For brighter fields, where the resid- 
ual also includes interstellar contributions, we assumed it also follows 
a k~ l power spectrum, but with a greater normalisation (i.e., cr s ) to take 
into account these effects. 



The FIRAS data provide a reference for the dust emission spec- 
trum of the diffuse ISM. The average FIRAS SED of the high- 
Galactic latitude sky used by Compiegne et al. (2011) to set the 
diffuse ISM dust properties in the DustEM model is shown in 
Fig. 17. Also shown in this figure are the current results (red 
symbols) for the average of the emissivities for the LVCs of our 
sample at 353, 545, 857, 3000, and 5000 GHz. The IRAS and 
Planck data points are found to be fully compatible with the dif- 
fuse ISM FIRAS spectrum, showing that our sample is represen- 
tative of the diffuse dust emission at high Galactic latitudes. We 
have fit model parameters to both data sets independently using 
a modified black body function: 



e v = Iy/Nm = Ko(v/voy/j,m u By(T), 



(9) 



where B V (T) is the Planck function, mu the mass of hydrogen, 
ji the mean molecular weight and kq is the opacity of the dust- 
gas mixture at some fiducial frequency vo. The higher-frequency 
(60 fim) IRAS datum is not used in the fit due to contamination 
by non-equilibrium emission from stochastically-heated smaller 
grains (Compiegne et al. 2011). The top panel of Fig. 17 shows 
the data divided by the model to better display the quality of the 
fit. The values found for T and J3 with the two data sets are in 
close accord. This analysis shows that the local ISM SED can be 
fit well with T = 17.9 K and J3 = 1.8. 

5.3. The spectral energy distribution 

Fig. 16 shows the dust SED for each field and Hi component 
separately, with the error bars computed using the Monte-Carlo 
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Fig. 16. SEDs from the emissivities of LVC (black) and IVC (blue) components for all the fields in our sample. For each H i 
component in each cloud, the solid line is the modified black body fit using 353, 545, 857, and 3000 GHz data. 



simulations. We first note that all the LVC and IVC SEDs are at a 
similar level showing that the power emitted per H is comparable 
between fields. Second, in many cases, the SED of IVCs peaks 
at a higher frequency than for the LVCs. That could be caused 
by a higher temperature or a larger abundance of smaller grains 
in the IVCs. 

Like for the FIRAS comparison, each SED was fit using a 
modified black body function (solid lines). It has been shown 
by several studies (Dupac et al. 2003; Desert et al. 2008; Shetty 
et al. 2009; Veneziani et al. 2010) how difficult it is to estimate 
separately T and /3 for such an SED fit. These two parameters are 
significantly degenerate; their estimate depends greatly on the 
accuracy of the determination of the error bars on the SED data 
points and on the correlation of the errors between frequencies. 
The estimate of T and j3 also depends on the spectral range used 
to make the fit; for a typical T = 1 8 K dust emission spectrum, 
the Rayleigh- Jeans range of the black body curve, where j3 can 
be well estimated, corresponds to v < 375 GHz (A > 800 jim). 

Because of the above caveats relating to simultaneous fits of 
T and J3, and because the FIRAS spectrum of the diffuse ISM 
is compatible with ft - 1.8 (in this case the large number of 
data points and the broad frequency coverage over the peak of 
the SED give more confidence in the value of J3 obtained), we 
have carried out SED fitting assuming a fixed ft = 1.8. This pro- 
vides a direct way to compare not only with the FIRAS spec- 
trum and the DustEM model, but also with similar analyses 
with Planck data on molecular clouds and in the Galactic plane 
(Planck Collaboration 201 lo,q,u) that used the same convention. 
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Fig. 17. Bottom panel: Black points show the FIRAS spectrum 
of the diffuse ISM (Compiegne et al. 2011). The red points are 
the average of the IRAS or Planck emissivities for the local com- 
ponents of all our fields; the uncertainty is the variance of the 
values divided by ^^N. The solid curves are modified black body 
fits to each spectrum - the 5000 GHz point was excluded from 
the fit as it is dominated by non-equilibrium dust emission. Top 
panel: Same as bottom panel but each data set is divided by its 
modified black body fit. 



5.4. Dust properties 

The modified black-body fit (Eq. 9) provides information on the 
properties, such as T, of the dust in each H i component. A useful 
quantity used below is the emission cross-section of the interstel- 
lar material per H: 

cr Q = K (v/v fidmn = r/Nuh (10) 

It is simply the prefactor to the Planck function in Eq. 9. In what 
follows we adopt vo = 1200 GHz (Ao = 250 /mi) to compare 
directly with the value of <x e (1200) = r/Nm at 250 yum given by 
Boulanger et al. (1996). 



A key quantity is the luminosity per H atom L (in W/H) emit- 
ted by dust grains (equal to the absorbed power) computed by 
integrating the SED over v: 

L = J 47iKo(y/yofj^m n B v (T)dy. (11) 

Complementing the actual SEDs in Fig. 16, Fig. 18 shows 
the derived values of <x e (1200), T and L plotted against the ve- 
locity of each H i component (LVC and IVC) in our sample. 

The average emission cross-section for the 14 LVC compo- 
nents of our sample is 1.0 ± 0.3 x 10" 25 cm 2 , in good agreement 



21 



Planck collaboration: Planck Early Results. XXIV. Dust in the diffuse interstellar medium and the Galactic halo 



2.0 
1.5 
1.0 

0.5 
0.0 



24 - 
22 
20 
18 



16 
0.6 

0.5 

0.4 

0.3 
0.2 




-80 



-60 



-40 -20 

v (km s~ 1 ) 



20 



Fig. 18. Values of cr e at 1200 GHz (250 //m), T and L versus the 
average velocity of each H i component. These dust parameters 
were estimated from the SED fit over the range 353 to 3000 GHz 
to a modified black body with j3 = 1.8. Black and blue sym- 
bols are for LVC and IVC, respectively. The different symbols 
in blue represent IVCs associated with specific complexes: IV 
Arch (dot), IV/LLIV Arch (triangle), Complex K (square) and 
PP Arch (star). In each panel the dotted line represent the value 
for the diffuse ISM obtained with the high-latitude FIRAS spec- 
trum (see Fig. 17). Error bars are given for each data point, some 
being smaller than the symbol size. 



with the value of 1 x 10 25 cm 2 obtained by Boulanger et al. 
(1996). The scatter of cr e (1200) for the LVCs (30%) is not the 
result of errors. 

The emission cross-section for the IVC components is dif- 
ferent, often 50% lower compared to the LVCs. There appear to 
be differences among the IVCs too, perhaps related to the fact 
that they belong to different IVC complexes. All fields in our 
sample overlap with the Intermediate Velocity (IV) Arch, ex- 
cept MC which is in the southern Galactic sky and belongs to 
the PP Arch and BOOTES which is part of Complex K (Kuntz 
& Danly 1996). The North Celestial Loop also overlaps spa- 
tially with the Low-Latitude Intermediate Velocity (LLIV) Arch 
(Kuntz & Danly 1996), an Hi feature at slightly less negative 

velocity ( 50 kms" 1 ) than the IV Arch ( 75 kms" 1 ). The 

fields UMA and UMAEAST contain clumps identified by Kuntz 
& Danly (1996) as being part of the LLIV Arch (specifically 
LLIV1, LLIV2 and LLIV3). The field POL also contains emis- 
sion that can be attributed to the LLIV Arch. These different 
complexes are identified with separate symbols in Fig. 18. The 



outliers are POL, UMAEAST (high) and MC (low). Excluding 
these cr e = 0.5 ± 0.2 x 10" 25 cm 2 for the rest. 

The values of T for the LVCs (T = 17.9 ± 0.9 K) are in 
accordance with that obtained from the FIRAS spectrum at high 
Galactic latitude (T = 17.9 + 0.1 K; Fig. 17) and also close to the 
17.5 K found by Boulanger et al. (1996) assuming p = 2. Like 
for <r e , a systematic difference in T is found between LVCs and 
IVCs, even though it is less statistically significant. On average 
the IVC group of clouds 6 has T = 20.0 ± 1.0 K, a value greater 
than in the local ISM at the 2.1 cr level. Note the different IVC 
complexes as well. 

Regarding L the striking result here is the small variation 
observed over all the fields and Hi components. Combined to- 
gether, the LVCs and IVCs have < L >= 3.4 ± 0.6 x 10" 31 W/H, 
representing a variation of only 20% over all clouds. The fact 
that L is rather constant over all fields and H i components can 
also be appreciated in Fig. 16, where all SEDs are at about the 
same level. The small variation of L is surprising as it indicates 
that the power absorbed by dust is also rather constant, possibly 
suggesting constancy of the radiation field across all fields even 
at the distance of IVCs. 

We have made the same analysis with a floating J3 to evalu- 
ate the robustness of the results shown in Fig. 18. The greatest 
impact of a floating ft is on the uncertainty of T, and to a lesser 
degree on cr e . The values of L are particularly insensitive because 
the modified black-body still serves effectively as an interpola- 
tion function between the measured data points. Even though the 
uncertainties on T and cr e increase with a floating J3, the trends 
seen in Fig. 18 are still observed, even more pronounced. We 
conclude that using ft = 1.8 is a reasonable and conservative 
approach. 



5.5. SED of the residual emission 

As seen in Figs. 4 to 11, the residual emission once the Hi 
model is subtracted from the Planck and IRAS maps exhibits 
significant spatial coherence which reproduces from frequency 
to frequency. In order to estimate the SED of this residual emis- 
sion in each field, we carried out a linear regression analysis be- 
tween the residual map at each frequency and the residual map 
at 857 GHz. The resulting slope of the regression (in units of 
MJy sr" 1 / MJy sr" 1 and equal to 1 at 857 GHz, by construction) 
can be used to estimate the shape of the spectrum of this residual 
emission. The results are presented in Fig. 19. 

As for the SEDs of the H i-correlated emission, the SEDs of 
the residual are well fit by a modified black body function. We 
note a significant difference in the SED shape between low H i 
column density fields and brighter regions. The SEDs of brighter 
fields, where the residual emission is likely to come from molec- 
ular gas (Sect. 6.1), could be fit with /3 = 1.8 but with a slightly 
lower temperature (T = 16.1 ± 0.6 K) than the LVC components. 
On the other hand, the SED of fainter fields, where the residual 
emission is most probably dominated by CIB anisotropies, could 
not be fit with p = 1.8. It is better described with T = 18.6 ± 0.9 
and ft = 1 . 1 ± 0. 1 . This should be taken only as a convenient fit- 
ting function, nothing physical. The significantly different value 
of p found here is probably the result of the complex composite 
nature of these fluctuations coming from the combination of all 
galaxies along the line of sight in slices of redshift which change 
with frequency (Planck Collaboration 201 In). 



excluding POL, UMAEAST and MC 
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Fig. 19. Correlation coefficient of the residual emission with the 857 GHz residual. Similar to Fig. 16, the solid line here is the 
normalised modified black body SED fit using 3000, 545 and 353 GHz data. We assumed a fixed J3 = 1.8, except for the six most 
diffuse fields (MC, G86, AG, Nl, SP and BOOTES) where the residual emission is dominated by the CIB fluctuations. In these 
fields J3 = 1.8, typical of diffuse Galactic dust emission, does not provide a good fit; these SEDs require a smaller value of f3 (see 
text). 



6. Discussion 

6.1. The Hi-H 2 transition 

One possible contribution to the residual emission is dust asso- 
ciated with ionized gas. The WIM has a vertical column density 
of about 1 x 10 20 cm" 2 (Reynolds 1989; Gaensler et al. 2008), a 
significant fraction of the total column density in the most dif- 
fuse areas of the sky. Planck Collaboration (201 In) showed that, 
once the emission correlated with H i is removed from the Planck 
353, 545 and 857 GHz data in faint fields, the residual emission 
has a power spectrum compatible with k~ l , much flatter than any 
interstellar emissions. We also showed that the amplitude of the 
residual in faint fields is constant from field to field (see Fig. 15), 
compatible with an isotropic extra-galactic emission. These are 
strong indications that dust emission associated with the WIM, 
and not correlated with H i, is small. 

In the eight fields with bright cirrus, the residuals to the IR- 
H i correlation are skewed toward positive values, larger than the 
amplitude of the CIB fluctuations. The most straightforward in- 
terpretation is that these positive residuals trace dust emission 
within H2 gas. This interpretation is reinforced by the fact that in 
the brightest fields (e.g., UMA, UMAEAST, and POL) we have 
checked that the residual emission is very well correlated with 
CO emission (Dame et al. 2001; Planck Collaboration 201 lo), 
confirming the previous study of Reach et al. (1998) in the North 
Celestial Loop region (see their Fig. 11). The lower dust tem- 
perature estimated from the SED of the residual emission in all 
these bright fields (see Fig. 19) is also reminiscent of what is 
observed in molecular clouds (Planck Collaboration 201 lu). 

In the following we assume that the submm excess emission 
provides a way to estimate the molecular gas column density 
within the local ISM (LVC) component 7 . In order to estimate 



7 This is not true for the DRACO field where the excess emission is 
likely to be dominated by Intermediate-Velocity gas which has accom- 
panying CO emission (Herbstmeier et al. 1993). For this reason this 
field was discarded in the current analysis. 



the fraction of H2 in our fields we carried out an analysis using 
the 7g57 dust maps and the following equation explained below: 



N 



LVC/ 



^857 _ 



IVC fylVC 
c 857 HI 



- 3 



HVC /yHVC 



857 



'857 



^LVC 
fc 857 



(12) 



To concentrate on the local (LVC) gas, we removed the IVC and 
HVC-correlated emission and the constant term from 7 8 5 7 , using 
the above results of the linear regression in each field. Assuming 
the dust emissivity is the same in the molecular gas as in the 
atomic gas from which it formed, we divided this by 6^ c for 
each field to produce estimated maps of the total column density 
A^ VC/ for the local/low- velocity gas. We note that this map still 
includes the fluctuations of the CIB (C857) that act like a noise 
on the true total column density A^j vc : 



aiLVC/ _ aiLVC , r 1 xvc 



(13) 



We of course have a map of the LVC atomic column den- 
sity (A^y c ) from the GBT measurements. Therefore, we are able 
to compute an estimate of the molecular column density map 
/yLvc = (#lvc/ _ n^ c )/2 and then calculate an estimate of 
the fraction of mass (or fraction of H nuclei) in molecular form, 
/(H 2 ) = 2N™ C /(2N™ C + N u/C) pixel by pixel xhe results 

obtained by combining all pixels in all fields, except DRACO, 
are plotted in Fig. 20. To produce this figure we binned the data 
in A^ VC/ and within each bin examined the PDF of /(H2), find- 
ing its median (dark symbol) and half-power points (error bars, 
which are not necessarily symmetrical). 

On this figure are also plotted estimates of /(H2) using FUSE 
data by Gillmon et al. (2006); Wakker (2006) (squares - high- 
latitude lines of sight) and Rachford et al. (2009) (triangles - 
Galactic Plane lines of sight). The combination of these two 
datasets covers the column density range probed in our analy- 
sis. 

The results plotted in Fig. 20 reveal an increase in /(H2) be- 
ginning at A^j vc ~ 3 x 10 20 cm~ 2 and reaching 0.8 for Nu ~ 
4 x 10 21 cm" 2 . Note that we are not sensitive to the much lower 
values of /(H2) found by FUSE at low column densities. For 
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1.0 




N H (10 20 cm- 2 ) 

Fig. 20. Fraction /(H2) of hydrogen that is in molecular form 
in local gas/LVC, calculated from emission excess relative to the 
linear correlation (see text), versus the total column density es- 
timated using the 857 GHz dust emission transformed into gas 
column density using the emissivities e v . Black points show the 
median value of /(H2) in bins of Nh computed using all lines of 
sight in our sample. The error bars show the half- width at half 
maximum of the PDF computed in each N H bin. Green symbols 
show the results obtained using UV absorption data from high- 
latitude surveys (Gillmon et al. 2006; Wakker 2006, squares) 
and toward O stars on lines of sight closer to the Galactic Plane 
(Rachford et al. 2002, 2009, triangles). 



pixels with A^j vc < 3 x 10 20 cm 2 the dispersion of our estimate 
of /(H2) is due to the fluctuations of the CIB. For these pixels, 
A^ vc « N^[ c and therefore N^ c « C 857 / 2 6^ 7 C . Because C 857 
has a mean of zero, it produces both positive and negative values 
of/(H 2 ). 

Note that the UV observations toward O stars in the Galactic 
Plane (triangles) give some /(H 2 ) values at the same level as we 
find, but also some much lower values for a given column den- 
sity. This suggests that the UV observations are somewhat af- 
fected by dumpiness and/or are sampling qualitatively different 
lines of sight than ours at high latitude. Indeed Wakker (2006) 
emphasized that it is not straightforward to relate the H 1-H2 tran- 
sition to physical properties of the interstellar gas, like density in 
the molecule-producing environment, or even a physical thresh- 
old in column density required for molecule formation, because 
of the summing over different environments along the line of 
sight. A corollary is that there can be quite different values of 
/(H2) for the same Nu, as is seen in the figure, and so compari- 
son of UV data with the complementary infrared/submm analy- 
sis is of great interest. 

The overall correspondance between the trends seen inde- 
pendently in the Planck and FUSE results supports our inter- 
pretation of the submm excess being caused by dust associated 
with molecular hydrogen. We would call this medium "dark gas" 
if it were not detected via CO (Planck Collaboration 201 lo). 
Although there are lines of sight where /(H 2 ) is quite high, 
summed over all lines of sight in our survey, our study shows 
that the excess emission at 857 GHz that is not correlated with 
Hi is only 10% of the total emission. Thus the fraction of the 
hydrogen gas mass that is in molecular form in this sample of 
the high latitude diffuse interstellar medium in the solar neigh- 
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Fig. 21. Temperature vs emission cross-section at 1200 GHz 
(250 yum) estimated from modified black-body fit with (3 = 1.8 
to data from 353 to 3000 GHz. Black is for LVCs, blue IVCs: IV 
Arch (dot), IV/LLIV Arch (triangle), Complex K (square) and 
PP Arch (star). The solid line respresents a constant emitted lu- 
minosity L corresponding to the diffuse ISM reference values 
(o- e = 1 x 10" 25 cm" 2 and T = 17.9 K - dotted lines). 



bourhood is quite low. Planck Collaboration (201 lo) estimated 
a value of 35% for the entire high-latitude sky above 15°, which 
includes higher column density lines of sight (this introduces yet 
another factor, H 1 self-absorption, to make the medium dark). Of 
this 35%, about half is traced by CO, leaving about 20% as "dark 
gas" not traced by CO or H 1. 

The structure and nature of the diffuse molecular gas can 
be studied using the maps of residual (excess) dust emission. In 
SPIDER (see Fig. 9), an intermediate column density field where 
there is very little CO emission detected (Barriault et al. 2010), 
we observe coherent filamentary structures in the residual map 
that we interpret as the presence of dust in diffuse H 2 gas without 
CO. We have checked that the structures cannot be accounted for 
by an underestimate of the 21 -cm line opacity correction, being 
present even with T spin as low as 40 K. In this field the submm 
residual provides a way to map the first steps of the formation of 
molecules in the diffuse ISM. 

6.2. Evolution of dust 

6.2.1 . Variations of the big grain emission cross-section 

Interstellar dust evolves through grain-grain and gas-grain inter- 
actions. Fragmentation and coagulation of dust grains are ex- 
pected to occur in the ISM, modifying not only the grain size 
distribution but also the grain structure. The data described here 
provide important evidence for dust processing in diffuse local 
clouds and IVCs. As we will elucidate, the evidence foreshad- 
owed in Fig. 18 can be seen in Figs 21 and 22. 

Fig. 21 shows the values of T and <x e (1200) already shown 
in Fig. 18 but here as a scatter plot. The solid line shows the 
expected <x e (1200) as a function of T for a constant emitted lu- 
minosity L (normalized to 3.8 x 10" 31 W/H, the average diffuse 
ISM values for T = 17.9 K and o- e (1200)=1.0 x 10" 25 cm" 2 ). 
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Fig. 22. Top: dust emission per Nhi at 857 GHz versus the 
3000/857 GHz (100/350/mi) ratio. Local (black), IVC (blue - 
dot is IV Arch, triangle IV/LLIV Arch, square Complex K), 
residual (green). Solid line is the DustEM model for the dif- 
fuse ISM (Compiegne et al. 201 1), with radiation field variations 
from G = 0.1 to G = 5. Dashed line is the same model but with 
a relative abundance of VSGs four times higher than the stan- 
dard diffuse ISM value. Dotted lines gives the local ISM fiducial 
values (G = 1.0). Typical uncertainties are shown for each Hi 
component. Bottom: 5000/857 GHz (60/350 yum) ratio as a func- 
tion of the 3000/857 GHz (100/350 /mi) ratio. 



The emission cross-section <x e reflects the efficiency of ther- 
mal dust emission per unit mass. Dust that emits more efficiently 
will have a lower equilibrium temperature, the trend seen. Note 
that the dust in very different environments has close to the 
same integrated emission (L) and therefore is absorbing about 
the same power (solid line). The comparison with the values 
found in the Taurus molecular cloud by Planck Collaboration 
(2011u) (T ~ 14.5 K, = 1.8, cr e ~ 2.0 x 10" 25 cm" 2 and 
L = 2.3 x 10" 31 W/H) suggests that the T - cr Q anti-correlation 
extends to colder dust and denser environments, although the 
typical power absorbed is lower because of shielding. 

Because cr e depends only on the dust properties and on the 
dust-to-gas ratio, the 30% variation observed in local ISM (LVC 
- see Fig. 18), where we expect little variations of the metallic- 
ity, is interpreted as a genuine variation of the dust properties 
from one cloud to another. An alternative explanation would be 
the presence of different amounts of H 2 gas spatially correlated 
with the Hi. In this case the dust opacity <x e (1200) would be 
overestimated due to an underestimate of the actual gas column 
density. Although this scenario cannot be excluded formally, the 
data shown here does not support an increase of <x e (1200) with 
A/hi through the Hi-H 2 transition, as one might expect in this 
case. For this reason we favor an interpretation of the variations 
of <x e (1200) related to modifications of the grain properties. 



Fig. 22 complements Fig. 21 by comparing directly mea- 
sured emissivities (e v ) with the DustEM model of the average 
diffuse high latitude emission. This comparison, independent of 
any modified black-body fit, also shows strong evidence for dust 
evolution. The top plot of Fig. 22 shows the 857 GHz emissiv- 
ity versus the 3000 GHz to 857 GHz ratio. This is compared to 
a simple DustEM model prediction for constant dust properties 
and a variation of the radiation field strength from G = 0.1 
to G = 5 (G = 1 being the fiducial value). The prediction 
of DustEM is an increase of both 6357 and 63000/^857 with in- 
creasing radiation field strength, and therefore T. The former 
increases with G simply because of the increase of B V (T). The 
ratio 63000/6857 increases with T as the peak of the black-body 
shifts towards higher frequencies. The data points do not fol- 
low this trend showing clearly that the variations in the SEDs 
found here in LVCs and IVCs cannot be explained by local vari- 
ations of the radiation field strength. The data are consistent with 
a decrease of the dust emission cross-section (~ 6357) with tem- 
perature (~ 63000/6857), the same trend seen in Fig. 21. An evo- 
lutionary model in which dust structure changes due to aggre- 
gation (and the reverse process, fragmentation) is qualitatively 
consistent with these results: grains with a fluffy structure will 
absorb about the same amount of optical and ultraviolet radia- 
tion per unit mass as more compact grains, but compared to these 
more homogeneous and spherical grains they are more emissive 
at submm wavelengths because of their more complex structure 
and therefore cool more efficiently (Stepnik et al. 2003). 

6.2.2. Dust shattering in Intermediate Velocity clouds? 

The bottom plot of Fig. 22 shows the 5000 GHz to 857 GHz 
emissivity ratio as a function of the 3000 GHz to 857 GHz ra- 
tio for all Hi components and the residual emission SEDs. The 
locus from standard ISM DustEM models shows an effect of in- 
creasing colours with increasing G. The bright cirrus clouds of 
the LVC components (black symbols at the left end of the plot) 
and the molecular residuals (green) have colours in good agree- 
ment with the standard DustEM model, consistent with some 
variation of G < 1 . 

Both ratios are higher for the IVCs and the trend is offset rel- 
ative to the standard locus. The qualitative interpretation is that 
the emission in these higher- frequency bands is more contam- 
inated by non-equilibrium emission from an increased relative 
abundance of very small grains (VSGs) compared to the larger 
grains (BGs) in thermal equilibrium. To make this more quanti- 
tative, the dashed curve passing closer to the IVC data shows the 
colours for a DustEM diffuse ISM model with a relative abun- 
dance of VSGs four times higher than the standard value. 

Because of their non-equilibrium emission, VSGs dominate 
the emission for v > 5000 GHz and, according to the standard 
DustEM model and for the standard radiation field (G = 1), con- 
tribute to 30% of the emission in the IRAS 60 fim band even with 
an abundance of only 1.6% of the total dust mass. By contrast, 
the BGs make up 90.7% of the total dust mass, 14.2% carbon 
rich and 76.5% silicates. Therefore an increase of a factor four 
of the mass in VSGs is easy to accommodate by shattering of 
a small fraction of the BG dust mass. The fact that the emis- 
sion cross-section of IVCs is a factor of two lower than for the 
local ISM can be due to a decrease of the mass in BGs (i.e., 
the dust-to-gas mass ratio) and/or to a decrease in the emission 
cross-section of a single grain, the latter affecting the equilibrium 
dust temperature. For the average temperature found in IVQ 
(T = 20.4 K), and considering a variation of the dust-to-gas 
mass ratio from 1 to 1/2 the local ISM value, the radiation field 
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£353 


^545 


^857 


£3000 


£5000 


AG 


< 0.010 


< 0.029 


< 0.055 


< 0.071 


< 0.024 


BOOTES 


< 0.079 


< 0.258 


< 0.519 


< 0.289 


< 0.106 


DRACO 


< 0.021 


< 0.075 


< 0.168 


< 0.194 


< 0.054 


G86 


< 0.035 


< 0.096 


< 0.232 


< 0.342 


< 0.087 


MC 


< 0.003 


< 0.009 


< 0.017 


< 0.027 


< 0.037 


Nl 


< 0.016 


< 0.047 


< 0.094 


< 0.054 


< 0.018 


NEP 


< 0.014 


< 0.040 


< 0.084 


< 0.085 


< 0.076 


SP 


< 0.015 


< 0.041 


< 0.068 


< 0.040 


< 0.015 


SPC 


< 0.012 


< 0.040 


< 0.090 


< 0.060 


< 0.030 


SPIDER 


< 0.017 


< 0.061 


< 0.144 


< 0.113 


< 0.054 


UMA 


< 0.007 


< 0.021 


< 0.056 


< 0.053 


< 0.024 


UMAEAST 


< 0.017 


< 0.054 


< 0.135 


< 0.105 


< 0.028 


LVC 
IVC 
HVC 


0.045 + 0.006 
0.023 + 0.015 
< 0.013 


0.18 + 0.03 
0.09 + 0.04 
< 0.040 


0.52 + 0.09 
0.27 + 0.09 
< 0.085 


0.55 + 0.10 
0.66 + 0.06 
< 0.077 


0.10 + 0.02 
0.20 + 0.03 
< 0.034 



Table 5. Upper limit for the HVC emissivities in each field. For the emissivities of the LVC and IVC components, weighted average 
and standard deviation of the 14 fields. Last row: average HVC upper limit, excluding BOOTES, G86 and SPIDER for which the 
H 1 column density of the HVC is less than 10 19 cm" 2 inside the mask. All values are in MJy sr _1 /10 20 cm -2 



strength needed to explain these results would be G = 1-2, not 
an unreasonable range for IVCs. It is thus difficult to conclude at 
this point on the exact nature of the variation of cr e in IVCs. On 
the other hand these results are clearly compatible with a modifi- 
cation of the BGs (total mass and/or grain properties) and a sig- 
nificant increase of the abundance of VSGs. The details would 
be linked to the specific dynamical and shock history of this in- 
terstellar matter that is part of the Galactic fountain (Shapiro & 
Field 1976; Houck & Bregman 1990). 

6.3. Dust in high-velocity clouds 

Dust emission in HVCs is expected to be weak for several rea- 
sons. First, their great distance from the Galactic disk (typi- 
cally several kpc - Wakker 2001; van Woerden & Wakker 2004) 
implies a faint radiation field, hence a low dust equilibrium 
temperature and a low luminosity (Wakker & Boulanger 1986; 
Trifalenkov 1993; Miville-Deschenes et al. 2005). But the main 
reason comes from gas phase abundance measurements obtained 
from absorption lines of oxygen, silicon and iron in the far-UV. 
Because oxygen is mostly in the gas phase and because its ion- 
ization potential is similar to hydrogen, the ratio O i/H 1 provides 
the most robust measure of metallicity in diffuse clouds. Unlike 
oxygen, iron and silicon are depleted into dust grains in the lo- 
cal ISM (Savage & Sembach 1996); low ratios of Si/O and Fe/O 
thus trace the presence of dust. In most HVCs the O i/H 1 ratio 
shows that the metallicity is significantly sub-solar (of the or- 
der 0.1-0.5 solar, see Lu et al. 1998; Wakker et al. 1999; Richter 
et al. 2001; Tripp et al. 2003; Sembach et al. 2004; van Woerden 
& Wakker 2004; Fox et al. 2006; Collins et al. 2007), indicative 
of their extragalactic origin. In addition the Si/O and Fe/O ra- 
tios in Complex C are infered to be close to solar (Richter et al. 
2001; Collins et al. 2007) indicating little silicon and iron locked 
in dust grains. Together the overall low abundance of metals and 
the fact that Si and Fe seem to be mostly in the gas phase imply 
a low dust to gas ratio in HVCs. 

Nevertheless, because of the unknown enrichment history of 
HVCs that might set different intrinsic abundances than in the 
solar sytem 8 , the difficulty to estimate the effect of ionization on 



8 O and Si are mostly produced in Type II supernovae, Fe is produced 
primarly in Type la supernovae. 



Sin and Fen to infer the abundances of Si and Fe 9 , and the un- 
kown nature of dust in HVCs 10 , the link between the metallicity, 
the depletions and the amount of dust remains indicative. One 
illustrative related example, albeit for an IVC, is the finding by 
Richter et al. (2001) of only mild depletion of Si and Fe for the 
IVC IV- Arch which nevertheless has dust emission at a compara- 
ble level to the local ISM. Thus the search for dust emission from 
HVCs is an interesting complementary way to put contraints on 
the physical conditions and nature of matter in the Galactic halo. 

After unsuccesful attemps to detect dust emission from 
HVCs using IRAS data (Wakker & Boulanger 1986; Fong et al. 
1987; Bates et al. 1988; Desert et al. 1988), Miville-Deschenes 
et al. (2005) claimed a detection in Complex C using a combi- 
nation of Spitzer, IRAS and GBT data. Compared to these pre- 
vious searches, the analysis presented here extend greatly the 
number of fields, the surface covered and the wavelength cov- 
erage. The fields used here overlap with four HVC complexes: 
Complex C (DRACO, G86, Nl, NEP, SP, SPIDER, BOOTES 11 ), 
Complex M (AG), Complex A (SPC, UMA, UMAEAST) and 
the Magellanic Stream (MC). HVC-correlated emission is not 
detected significantly in any of these fields. 

At each frequency the emissivities found are compatible with 
zero within 3<x, with about as many negative as positive values. 
At first glance, examination of the emissivities in Table 2 shows 
systematic behaviour across frequencies that could be taken as 
encouraging evidence for an HVC detection in some fields, like 
AG. However, this does not necessarily bolster the statistical sig- 
nificance of the detection in a single pass band. This is because 
the dust emission is intrinsically highly correlated frequency to 
frequency, as is the CIB, and so given the generally small noise, 
the fits to the correlations with (the same) Nui are to first or- 
der just scaled versions of each other. This analysis establishes 
the main difficulties hampering the detection of dust emission 

9 the true abundance of Si and Fe are less than the ones infered from 
Si 11/H 1 and Fe 11/H 1. To take this bias into account, an ionization cor- 
rection is applied (Collins et al. 2003) that depends on gas density and 
radiation field strength, both not well constrained. 

10 the amount of carbonaceous dust is not constrained by gas phase 
abundances. 

11 BOOTES is on the side of Complex C. It has faint emission from 
-80 to -100 km s _1 that could also be associated with the IVC Complex 
K. 
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from HVCs, the contaminating foreground emission of LVCs 
and IVCs and the background CIB fluctuations which greatly 
increase the uncertainties in e^ vc (Sect. 4.5 and Table 4). 

In Table 5 we focus on the upper limits to the HVC emissivi- 
ties. They were computed using a Bayesian approach in the case 
of Gaussian errors with a bounded physical region (Feldman & 
Cousins 1998); we assumed that the true emissivities can only be 
zero or positive. Table 5 gives the value of e y for which there is 
99.87% probability (equivalent to a 3cr upper limit) that the true 
emissivity is < e y . These upper limits are compared in Table 5 
to the weighted average of the measured emissivites for the LVC 
and IVC components; any dust in HVCs has an emissivity, on 
average, no more than 15% of the value found in the local ISM 
at 857 and 3000 GHz. 

We considered the null hypothesis that the HVC emissiv- 
ity is zero (or a small uniform value) in all fields and find that, 
at any given frequency, the weighted standard deviation of the 
measured HVC emissivities (<x(e)) is somewhat larger than ex- 
pected given the uncertainties on each measurement. This is true 
even when considering only AG, Nl, MC and SP, the fields with 
significant HVC H i column density and almost no masking. For 
example, the probability of obtaining the measured cr(e) is only 
0.2% at 857 GHz for these four fields. One explanation could be 
that we have underestimated the uncertainties on each measure- 
ment; to increase the probability of the measured cr(e) to 10% at 
857 GHz, the uncertainties would need to be multiplied by 1.6. 
But even such a moderate increase of the uncertainties seems 
unlikely in the most diffuse fields considered here because the 
noise is dominated by the strongly constrained and well modeled 
CIB fluctuations (Planck Collaboration 201 In). A second expla- 
nation is that the cr(e) reflects a variation of the true e from field 
to field. For these four fields we verified that, if one allows for 
a variation of e from field to field following a uniform distribu- 
tion in the range between zero and the upper limit of each field, 
the probability of obtaining the measured cr(e) reaches 15% at 
857 GHz without having to increase the uncertainties on the in- 
dividual measurements. If this is the right explanation, the mea- 
sured value of cr(e) implies that e is small but non-zero in some 
fields. 

Finally, there have been many interesting discussions of the 
possible relationships between IVCs and HVCs (e.g., Albert & 
Danly 2004; Wakker 2001). Here the relative emissivity is of 
some relevance; as our data show, H i-correlated emission is 
readily detected in IVCs, but not in HVCs. A specific example 
is provided by our AG field which encompasses the cloud MI in 
HVC Complex M (Wakker 2001). The low HVC emissivity in 
contrast to the fairly normal IVC emissivity (Table 2) does not 
support any direct relationship between MI and the intermediate 
velocity gas in the IV Arch in this direction (Wakker 2001). 



7. Conclusions 

We have presented results of a first comparison of Planck and 
IRAS with new Hi 21 -cm line GBT data for 14 high Galactic 
latitudes fields, covering about 825 deg 2 on the sky. Using the 
velocity information of the 21 -cm data we made column den- 
sity maps of the LVC, IVC, and HVC in each field. By correlat- 
ing these H i maps with the submm/infrared dust emission maps 
we estimated the distinct dust emissivities for these three high- 
latitude components and made corresponding SEDs from 353 to 
5000 GHz. 

On average, the dust SED of the local ISM (LVC) deduced 
from the IRAS and Planck data is compatible with that from the 



FIRAS data over the high-latitude sky and is well fit with a mod- 
ified black body with parameters ft = 1.8 and T - 17.9 K. On 
the other hand, even though the energy emitted by dust is rather 
constant in our sample, we report significant variations of the 
dust SED shape, compatible with variations of the dust temper- 
ature anti-correlated with the emission cross-section. We inter- 
pret these variations as a signature of active evolution of the dust 
grain structure through coagulation and fragmentation in the dif- 
fuse ISM. 

For faint cirrus fields with average H i column density lower 
than 2 x 10 20 cm -2 , the residual FIR-submm emission, after re- 
moval of the H i-correlated contributions, is normally distributed 
with a standard deviation compatible with the expected level of 
CIB fluctuations. The SED of the residual is unlike any inter- 
stellar component. For such diffuse fields we also show that the 
interstellar dust emission is dominated by the contribution from 
atomic gas. 

For fields with larger Hi column density there are signifi- 
cant FIR-submm emission excesses that, for the brightest fields, 
follow the structure of the CO emission. For intermediate col- 
umn density regions, the residual emission shows coherent spa- 
tial structures not seen in CO, revealing the presence of H2 gas. 
The SED of this residual shows that it is slightly colder than the 
dust in the local H 1. There is a lot of variation in /(H2), the frac- 
tional mass in the form of H2, but over the whole survey of these 
diffuse fields the fraction is only 10%. 

We report strong detection of dust emission from all IVCs 
in our sample. The dust emission cross- section is typically two 
times lower than for the local ISM and the relative abundance 
of small grains having non-equilibrium emission is about four 
times higher than normal. This evolution of the dust properties 
is indicative of dust shattering in halo gas with the dynamics of a 
Galactic fountain. We also find that, compared to the local ISM, 
several of these clouds have a higher dust temperature which 
could also be the result of globally smaller grains. The total en- 
ergy emitted by dust in IVCs is comparable to what is observed 
in the local ISM suggesting similar radiation field strength. 

Finally, we have attempted to detect HVC-correlated dust 
emission. We show why this is very challenging, because of the 
uncertainties induced by foreground contamination and the CIB 
anisotropies. The average of the 99.9% confidence upper lim- 
its on the emissivity, 0.15 times the average emissivity found in 
the local ISM at 857 and 3000 GHz, is compatible with the low 
level of dust emission expected from HVCs based on their low 
metallicity and relatively low depletion of elements like Si and 
Fe. 
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Appendix A: Noise of the H i components 

The analysis presented here relies on a correlation between far-infrared/submm 
data and the H i components deduced from 21 -cm emission. In order to estimate 
properly the uncertainties on the correlation coefficients the noise on the H i 
components needs to be evaluated. We have made this estimate using two 
different methods which address different contributions to the noise (see 
Boothroyd et al. 2011). 

First we applied the traditional method in which the end channels of the 21 -cm 
spectrum, where there is no emission. The velocity channels are independent, 
i.e., the spectrometer resolution is better than the 0.8 km s _1 channel width. A 
large number of end channels can be used to compute a map of the standard 
deviation of the noise at each (x, y) position: 



^ 2 (x,v,v)-i-(f T b (x,y,v) 



(A.l) 
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where N§ = v2 - vl + 1 is of the order of 200 channels. For the GBT data, this 
is a very flat map, whose typical value is the same as the standard deviation of a 
single end channel. This is about 0.17 K for a single visit to a field. Some fields 
were mapped two or three times. 

The uncertainty map SNm(x, y) for the column density of a given H i 
component summed over N channels is 



Field 



5Nm(x,y) [cm 



: A VNS T b(x,y)Av, 



(A.2) 



where Av is the channel width in km s _1 . For N ~ 60, this amounts to 
dNni ~ 0.02 x 10 20 cm -2 for a single visit. This assumes, as is usual, that the 
noise properties estimated using the end channels are representative of the noise 
in channels summed to build the column density maps. Actually, the channel 
noise scales as roughly (1 + T^/20K) and this can be taken into account. 
We developed a second method which exploits the fact that the H i 
brightness-temperature cubes were built by averaging independent data taken in 
different polarisations (called XX and YY). Each polarisation observation can be 
reduced separately, which requires separate baseline estimates for the spectra in 
the XX and YY cubes. Baseline fitting is another source of error. Taking the 
difference between the XX and YY cubes removes the common unpolarized 
H iemission, leaving in each velocity channel only the thermal noise and 
systematic effects from the baseline subtraction. The difference cube is 



A(x,y,v) = (T X x(x,y,v) - T Y y(x,y, v))/2, 



(A.3) 



where we divide by two to give the same statistics as in the average of the cubes. 
We make a map of the column density differences (divided by two) over the 
same channel range as for Nui- 
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Table B.l. Planck and IRAS noise levels (in MJy sr" 1 ) for each 
field and at each frequency. The value given here is the noise 
level in the map that has been convolved to the 9.4' GBT resolu- 
tion. 



AN m (x, y) = A Av ^ A(x,y, v). 



(A.4) 



The estimate of SNui is then the standard deviation of this map. 
This method gives an uncertainty typically 1.3 times the simple estimate using 
end channels (it can reach up to 2 times for the brightest components) because it 
includes baseline uncertainties and the increase of thermal noise with signal. 
The column density uncertainties evaluated with this second method for each 
component in each field are given in Table 1). There are additional errors from 
the stray radiation correction that can affect the local and IVC components. The 
total uncertainty including this additional contribution can be estimated in a 
similar way to the second method if there are two or more visits to compare 
(Boothroyd et al. 201 1). Because in cases with multiple visits these are found to 
be not much larger, and because the uncertainties in A/hi are not the major 
source of uncertainty in the correlation analysis (Sect. 4.5), adopting the values 
in Table 1 is satisfactory for our purposes. 



Appendix B: Estimating the noise of the Planck and 
IRAS maps 

To estimate the noise level of the IRAS and Planck maps we use the difference 
between independent observations of the same sky. For IRAS, the original ISSA 
plates (Wheelock et al. 1993) were delivered with three independent set of 
observations (called HCONs), each obtained at different periods during the life 
of the satellite. Each HCON has its own coverage map. The ISSA final product 
is the coverage- weighted co-addition of the three HCONs. The IRIS product 
(Miville-Deschenes & Lagache 2005) also contains the three HCONs. 
A similar approach was used to estimate the noise of each Planck map. For each 
pointing period Planck scans the sky about 50 times. The Data Processing 
Centre delivered two maps made with data from the first half and second half of 
each pointing period, respectively. In this case the number of samples used to 
estimate the signal at a given sky position is the same in the two maps (i.e., the 
coverage map is identical for the two maps). Again we used the difference 
between these two maps of the same region of the sky to estimate the noise. 
Here we want to estimate the average noise level SI V of an IRAS or Planck map 
I v used for the analysis given a difference map 



A v (x,y) = (I v , 1 (x,y)-I v , 2 (x,y))/2 



(B.l) 



obtained from independent observations, each I v j map having its own coverage 
map N V j. In the general case (Miville-Deschenes & Lagache 2005), 8I V is 
obtained by taking the standard deviation of the map 



A v (x,y) = A v (x,y) 



4N v , l (x,y)Ny, 2 (x,y) 
Ny(x,y)(Ny^x,y) + N va (x,yy 



(B.2) 



In the case of Planck, N Vt \ = N Vf 2 = N v /2 and so A' v = A v and the standard 
deviation of the average map is the same for the difference map, as expected. 



Because the ISSA plates are the combination of three observations with 
different coverages such a simplification is not possible and the above more 
general equation has to be used. 

In our analysis we used the IRAS and Planck maps convolved to the 9.4' 
resolution of the GBT. Therefore, the appropriate 6I V is obtained by first 
convolving the A v (x,y) map and then taking the standard deviation. These are 
the noise levels tabulated for each field in Table B.l. 
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